VARIATIONAL EQUATIONS OF LAGRANGIAN SYSTEMS AND 
HAMILTON'S PRINCIPLE 

H. N. NUNEZ- YEPEZW, JOAQUIN DELGADCK 2 ), AND A. L. SALAS-BRITO< 3 ) 

f - *^ ■ Abstract. We discuss a recently proposed variational principle for deriving 

r™ "\ ' the variational equations associated to any Lagrangian system. The principle 

^vj t gives simultaneously the Lagrange and the variational equations of the system. 

We define a new Lagrangian in an extended configuration space — which we 
call DAlambert's — comprising both the original coordinates and the com- 
patible "virtual displacements" joining two solutions of the original system. 
The variational principle is Hamilton's with the new Lagrangian. We use this 
formulation to obtain constants of motion in the Jacobi equations of any La- 
grangian system with symmetries. These constants are related to constants 
in the original system and so with symmetries of the original Lagrangian. We 
cast our approach in an intrinsic coordinate free formulation. Our results can 
be of interest for reducing the dimensions of the equations that characterize 
perturbations in a Lagrangian control system. 
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1. Introduction 



i-£h ■ Extremum principles have a great importance in dynamics and in control [0,0, 

'■jp H j J3|. They can be used to attain general formulations of the theory and can be 

^£1 • of assistance in the solution of particular problems p|,||. The purpose of this 

contribution is to discuss a recently proposed extremum principle which generates 
both the equations of motion and the so-called variational equations describing 
deviations between any two solutions of a Lagrangian system B . The variational 
principle is Hamilton's but with a different Lagrangian. Since we are dealing with 

J^> , Lagrangian dynamical systems, in which such "least action" principle is assumed to 

5^ ■ be valid, the variational equations describe deviations between extremal curves in a 

certain configuration space so they, slightly abusing the language, may be called the 
Jacobi variational equations of the Lagrangian system as we have done in g] . Please 
note that though we use a notation with a strong flavour of classical mechanics 
|7|, the formulation is not limited in any way and can be applied to non-linear 
evolution equations, like the Einstein, the Korteweg-DeVries, or the Kamdomtsev- 
Petviashvili equations ||— |Q. The intrinsic formulation of section 4 should make 
clear this matter. The variational equations can also be used to study the stability 
of solutions and as starting points to evaluate their Liapunov spectrum as we have 
succintly discussed in Q. They are, besides, describable in the vertical extension 
of the Lagrangian formalism jl2| . We should mention previous related work along 
the same lines 1l| and an older contribution with different purposes (14]. 
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The article is organized as follows: In section 2 we define the new Lagrangian, 
use it to formulate the variational principle whose extremals are solutions of the 
Lagrange and the variational equations of the system. In this section we also dis- 
cuss basic features of the formulation. In section 3, we address the existence of 
constants of motion in the variational equations and its relation with constants in 
the original system. In section 4 we cast our approach in an intrinsic language, 
making connection with the symplectic manifold that can be associated with the 
prolongation of the original Lagrangian. If you are interested mainly in the math- 
ematical setting of the formulation this is the section to read. In section 5 we give 
various examples in which the approach may be used. 

2. The variational principle 

Let us consider then a iV-degree of freedom dynamical system endowed with a 
Lagrangian function L(q a q a , t), a — 1 ... N; defined in the tangent bundle TQ of its 
configuration space Q, i.e. L : TQ —> M. From the Lagrangian we usually construct 



the action @,|| 



1.2 



(1) %(*)]= / L{q a q a )dt 

where we are using square brackets to indicate that S is a functional of the paths 
q a (t) joining two given points, q a {t\) and qafa) (a = 1, . . . , N), in Q at two fixed 
instants of time t\ and ti. The extremalization of S directly leads to the Lagrangian 
equations of motion of the system [equations (ph below]. Let us now consider a 
different configuration space D, comprising both the original configurations of the 
system plus all the possible "virtual displacements" joining any two extremal paths 
of (pi). With the help of L, we can define a new (the mathematical meaning of all 
this should become clear in sections 4.2 and 4.3) Lagrangian 7(q, q, e, c, t) f| as 

,„. , . . . dL . dL 

(2) ~f{q,q,e,e,t) = T —e a + -—e a 

oq a oq a 

here, as in all of the paper, the summation convention is implied for repeated indices. 
The function 7 : TD — » R, where TD is the tangent bundle of D, is central in our 
formulation. The 2./V-dimensional D'Alambert configuration space D is assumed, 
mostly in sections 2 and 3, to be coordinatized by (q a , e a ), a = 1, . . . ,N. The ./V- 
component object e — (ei, 62, • • • , ejv) stands for the displacement from an extremal 
path of (Q) to another, and k — (ei, £2, . . . , ejv) for its corresponding velocity. This 
means that e plays the role of the variational field associated with trajectories of 
the original system: q' = q + c and q' = q + c, where both q and q' are extremal 
paths 0], Q, H of the original action (|l|). A useful property of 7(q, q, c, e, t) is 
that it is an explicit function of time only when L is non-autonomous. 
The function 7 is used to define the new ( "displaced" ) action functional 

(3) E[q(t),e(t)]= I 7(q, q,e,e,t)di 

of the paths joining two given configurations (q\, e\) and(^2, £2) of the varied system 
between two fixed instants of time t\ and £2- We are calling "time" the parameter 
appearing in (||), but it can be any other useful parameter, for example, the arc 
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length, or even may be any finite set of parameters. The statement of the variational 
principle is just Hamilton's, that is, E is extremal when the system follows its actual 
path in DM 



(4) 



*£[q(t),e(t)]=0. 



The extremalization is done varying the path but mantaining the cndpoints and 
the time fixed. The conditions for the functional S[q(i),e(t)] to be an extremum 
are the 2N Euler-Lagrange equations [| , [Q , [jlij 



,,N; 



dt\de a ) de a ' dt\dq a ) dq a 

or, using the definition (|2|) in the preceding equations, we obtain the Lagrange 
equations of the original system 



(6) 



d_ f dL_ \ _ dL 
dt \dq a ) 8 

plus the linear variational equations: 



= 0, a=l,...,N, 



(7) 



M a bh + C ab e b + K ab € b = 0, 0=1,. 



,N. 



Equations (|7J) describe the evolution of the deviation, e, of a varied trajectory from 
an unperturbed one. The N x N matrices M, C and K, are 



(8) 



(9) 



(10) 



M ah = 



Cab — 



K ab = 



d 2 L 



s dq a dq b J ' 






d f d 2 L ^ 

dt \ dq a dq b J 


d 2 L 
dq a dqb 


d 2 L 

dqbdq a 


d f d 2 L \ 


d 2 L ' 




dt V dq a dq b J 


dq a dqb _ 


, a, b = 



,N. 



It shoud be clear that equations (|7]) are the variational equations of the original 
system [|j , || , g . One of the interesting properties of these equations is that the 
matrices M, C, and K, and, hence, the variational equations (pi), may be regarded 
as time-independent as long as L is autonomous. This property, that does not 
happen in the standard description, may have important consequences in analysing 
perturbations to optimal control problems || |T|] |T7J]. At this point, it is worthwhile 
to emphasize four important features: 

1. The function 7 : TD — > M. plays the role of a new Lagrangian describing the 
original system plus its response to "virtual displacements" or perturbations. 
It therefore is also useful for studying stability Q. 

2. Equations (o) and (m) are invariant under arbitrary point transformations (i.e. 
changes of coordinates in Q) : Q a — f a (q, t), a = 1, . . . , N, where the f a are 
functions, f a : D — > D, assumed to be invertible (i.e. det(df(q)) ^ 0, q E D) 
and at least C . As happens in any Lagrangian description, this result is 
proved by the very existence of the variational principle (||); p|,fl,|l7||. 
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3. The e-derivatives of 7 reduce to corresponding (/-derivatives of L 

dj dL d-y dL 

oe a dq a de a oq a 

4. As a consequence of the previous property, 7 is a first-order homogeneous 
function of the virtual displacements e a and velocities i a 

(12) 7 = ^T~ £a + -TT^ta- 

de a de a 

3. Symmetries and constants of motion 

One of the most important consequences of any Lagrangian description using 
a variational principle, is the close association between the symmetries of the La- 
grangian and the existence of constants of motion | J18| , that is, functions that are 
conserved along integral curves of the Lagrangian vector field. There is a control- 
theory version of such result, saying that every symmetry of a control system gives 
rise to a conservation law along biextremals H . This theorem makes the following 
results important for perturbed control systems. 

The association between symmetries and conservation laws is shared by the the 
function 7. The precise sense is the following: Any symmetry of L can be carried 
over to a symmetry of 7, where they imply the existence of constants of motion 
in the variational equations. Although a general proof can be given, we prefer to 
illustrate this result in what follows — though our discussion of what we have called 
inherited constants comes close to being an informal proof. 

3.1. Invariance under time translations. If the Lagrangian of the system is 
autonomous, then there exist a well-known constant of motion |fl] , H , 

(13) H=^q b -L. 

oqb 

H can, in certain instances, be identified with the energy of the system. 

In such autonomous case, 7 is autonomous too and the variational equations 
thus admit an analogous constant of motion — which can be derived in an strictly 
similar fashion to H JO. [[1,0] namely 

(14) h=^Lq b + ^e b -l- 

oqb oe b 

Using definition (||), the constant h can be recasted as || 

(?7 c?7 

(15) h = —q b - — e b , 

oq b oe b 



or, using (|13[), as 

, 9H. dH 

(16) h= ^—e b + ——e b . 

dq b dq b 

If we interpret H as the energy, a related interpretation for h is that it is the 
first-order energy change in going from a solution to a displaced nearby one. 
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3.2. Invariance under space translations. The invariance of a system under 
space translations, the property frequently called by physicists the homogeneity 
of space, usually manifest itself in the independence of the Lagrangian on certain 
coordinates. When the Lagrangian does not depend on a specific coordinate q s , the 
coordinate is said to be ignorable (respect to the original Lagrangian L), then its 
conjugate momentum is conserved. That is, if 

.._. dL d"f dL 

(17) — — = 0, then p s = —- — 777- is a constant. 
dq a de s dq s 

If q s is ignorable in L, then, as shown in equation (|17|), it is also ignorable in 7, 
and e s is ignorable too, so the momentum p s — conjugated to e a in 7 — is conserved 
(as it should be since it is conserved in the original system!). Furthermore, q s is 
also ignorable in 7, hence the quantity (the momentum, ir s , conjugated to q s in 7) 

, 10 , 9 7 d 2 L . d 2 L 

(18) n s = — - = eb + ». » €b, 

dq s oq s dqb dq s dqb 

is also a constant of motion in the variational equations. 

After having proved the existence of the constants h (equation pL4j ) and it (equa- 



tion 18), it is worth formulating a more general result. As we show next such 
theorem follows from the close relationship between solutions to equations (g) and 
constants of motion in the variational equations (f7|). We call such first integrals, 
inherited constants of motion. 

3.3. Inherited constants of motion. Let us consider any constant of motion, 
J(q, q), of the original set of equations (6). If we evaluate it on two nearby solutions 
of (6) , q and q' = q + e, separated by the Jacobi field c, the difference, j (c, e) = 
^(q'i <i') — ^(q ; q)i between these constant quantities is also trivially a constant, 

< 19 » *!Mi J <o'«- J <'»'«i= ' 

The constant, j[e], can be expressed as 

(20) j[ e ]=( ea+ 6c 

\ oq a dq a 

This is precisely the form of h in equation (Qq) . 

Equation ( |20| ) tells us how, given both a solution of equations (||) and any one of 
its constants of motion, we can obtain a constant of motion in equations (H) . Equa- 
tion (EG) can be directly proved to be a constant by computing its time derivative. 
This result establishes a direct relationship of constants in the variational equa- 
tions, like j, to constants in the original Lagrangian. Related results are discussed 
in §,ll|. 

Furthermore, as in some non-linear evolution equations the constants of motion, 
J[q(t)], are functionals (and not functions) of the solutions of (||), we need to 
pinpoint that the constant j[e(£)] becomes a functional of the Jacobi fields e(t), 
which must be given thus by 

(21) ^ = /w £(T)dT ' 
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where 5J[q(t)]/5q(r) stands for the functional derivative of the constant functional 

J[q(t)} gg. 

Having stablished the existence of the inherited constants (]2^) [or (El])], it should 
be clear that any symmetry of the original Lagrangian is reflected in the existence 
of another constant in the Jacobi variational equations. Thus, Noether's theorem 
also holds for the variational equations and can be used to reduce the dimensions 
of the variational system H;||. Related results are discussed in J7|, section 111]; 
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4. Intrinsic formulation 

In this section we formulate in an intrinsic manner our variational approach. 
First we consider an arbitrary vector field Y (at least C 1 ) defined on a manifold 
M and recall the constructions of the variational vector field T(Y) defined on the 
tangent bundle TM and the adjoint variational vector field T*(Y) defined on the 

cotangent bundle T*M. As it is well known @, Jl| , H) , || , || , T *^ YX ^ P osseses 
a natural structure: it is a Hamiltonian vector field with respect to the the canonical 
symplectic form on T*M. On the other hand the variational vector field T(Y) does 
not have a natural Lagrangian structure due to the fact that the Hamiltonian for 
T*(Y) is linear in the momenta and so there is no natural Legendre transform 
TM -> T*M defined. 

Our main interest is in the case of M = TQ, i.e. the tangent bundle of a 
configuration space, and Y = Yt, a Lagrangian vector field on TQ for the function 
L : TQ — > K. We recall that if L is not degenerate, the pull back of the canonical 
symplectic form uio in T*Q under the Legendre transformation C(L) : TQ — > T*Q, 
ujl = £*wo, makes (TQ,lul) into a symplectic manifold and the Lagrangian vector 
field Yl is Hamiltonian for the energy function H : TQ — > R. Also Yj, is a second 
order differential equation: tqJYl =1[Ot where tq : TQ — > Q is the projection. 

From what was said in the first paragraph of this section, the adjoint varia- 
tional vector field T*(Yl), defined on T*(TQ), has a natural Hamiltonian form. 
We want to dilucidate what the structure of the variational vector field T(Yl), 
defined on T(TQ), is. A naive guess is that T(Yl) should be Lagrangian and the 
Lagrangian should be the prolongation Jl|,|| L: T(TQ) -> R of L: TQ -> M: 
L(£) = d£„.(£)(£), for £ G T(TQ). A moment of caveat shows that this cannot be 
true. In fact, taking coordinates (ej», q, e, e) for T(TQ), and (q,q,p,ir) in T*(TQ), 
the prolongation is given by 

,_,_. dL . dL dL 

(22) —{q,q,e,e) = —(q,q)e+—{q,q)e 



aim ils ijegenuie iransiorm l 
(23) 


[1 K4) > 1 [1 K4) 

dL dL 
de dq 


(24) 


dL dL 
de dq 


but then the pull back of 9 - 


= pdq + 7r dq is 
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dL dL 

(25) 6 L = C(L)*9 = (poL)dq + (iroL)dq = —dq+—dq 

that is 9 1 = dL o ttq, there ttq : T(TQ) — > TQ is the projection. This in turn 
implies d9^ = (!). Thence the interpretation is untenable. 

4.1. The variational and adjoint variational vector fields. For setting straight 
the formalism, let us recall briefly the appropriate construction: Let M be a man- 
ifold and let t m ■ TM -> M, t* m : T*M -> M be the tangent and cotangent 
bundles, respectively. Given a vector field on M, that is a section Y : M — > TAf, 
the variational vector field associated to Y, herein denoted by T(Y), is defined 
as follows: Given m € M there exists a neighborhood B of m and e > such 
that the local flow $ t : B 4 (m) -> M is defined for t 6 (-e,e). The lift T($ t ) : 
TM| T -i(B) : 7r _1 (i3) — > TM then defines a local flow and its infinitesimal generator 
is the variational vector field T(Y). Similarly, the lifting to the cotangent bun- 
dle T*(4> t ) : TM\ p -i {B) : ir~ l {B) -> T*M defines a local flow and its infinitesimal 
generator is called the adjoint variational vector field T*(Y). 

The names are well suited since in local coordinates, (x, v) for TM and (x,p) in 
T*M (we think x, v as column vectors; p a row vector), the variational and adjoint 
variational vector fields are given respectively as 



x = f(x), x = f(x), 

v = Df{x)v; (26. a) p = -pDf(x) (26.6) 

Actually the adjoint variational vector field J2q.b) admits the Hamiltonian form 

where 

(29) H(x,p)=pf(x). 

This form is compatible with the standard symplectic form on T*M. As we men- 
tioned earlier in this section, there is no natural Lagrangian structure for (p6|.a) 

r — 1 / tt\ 

since the Legendre transformation (x,p) ► (x,v) is not defined due to the 

linear dependence of H on the momenta (|29j) . 

4.2. Prolongation of a system. There is a related concept of extending a given 
system known as its prolongation H , [ [15[ , p6| . The motivation comes from the fol- 
lowing argument: if the system is given by x — f(x), then by taking time derivatives 
along solutions one gets 

(30) 

x = Df(x)x = Df(x)f(x), x - D 2 f(x) ■ (Df(x)J(x) + Df(x) ■ (Df(x) ■ f(x)), etc. 
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Considering the first prolongation x = Df(x)f(x), it can be written as the first 
order system 

(31) v = Df(x)-f(x). 

A natural question is if this system admits a Lagrangian structure. 

Proposition 1. Consider the system x — f(x) and suppose that -J^- = -J^- , that is 
Df(x) be symmetric. Then the first prolongation of the system ( pip is Lagrangian 
for the function 

(32) L{x,v) = \\v\* + \\f{x)\. 

The proof is elementary and will be omited. □ 

The symmetry condition can be rephrased in terms of the graph 

(33) G f = { (x, f(x) T ) | x G Dom /} C R n X (R n )* 

as being Lagrangian p4| . ]l5[ | (here the notation f(x) T means simply that f(x) is 
viewed as a covector). 

If /: M -> N, its prolongation ||Jl| is the map /: T(M) -> T*(N) defined as 
follows: For £ 6 T X M and a £ T1,JM), let ( , ) denote the natural pairing, then 

(34) /(O(o) = (/*x(0.a>- 

In local coordinates (&,«) £ I" x R™ for TA/J and (y,p) £ R" x (R")* for T*(iV) 
the first prolongation is given by 

(35) y = /(s), 

(36) P = v T Df(xf. 

In the case of a function /: M — > R, TJ (:r sR ~ R and the prolongation can be 
identified with the differential of the map df : TM — » R 

(37) /(&) = 4f x (&). 

4.3. Variational equations of Lagrangian systems. We now consider the case 
of M — TQ for some configuration space Q, and L: TQ — > R a non degenerate 
Lagrangian for the vector field Y on TQ. The prolongation is then the a map 
i: T(TQ) —¥ R. As was mentioned at the beginning of this section, the varia- 
tional vector field T(Y) cannot be viewed as the Lagrangian vector field for the 
prolongation L in an obvious way since the pullback of the canonical 1-form in 
T*{TQ) to T(TQ) by the Legendre transform C{L) of the prolongation L is an 
exact differential. 

Instead we consider the D'Alambert phase space consisting of pairs (q,e) where 
q is a given configuration and e is a "virtual displacement" . More formally we 
consider a subbundle , D C Q then TD an be embeded in T(TQ) under the map 
(q, e, (j, e) — ► (q, q, e, e). Now view the prolongation as a map 7: TD — > R, that is 
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dL dL 

(38) 'y(q,e,q,e) = —(q,q)e+—(q,q)e 

(in fact TD is the natural domain of definition of the prolongation, since the "virtual 
displacements" are defined by choosing the subbundle) . 

Notice that the embeding a is well defined since the e's are elements of the 
tangent space T q Q, and that D is not necessarily a tangent bundle, that is, it is 
not necessarily an integrable distribution. Thus we can consider even holonomic or 
nonholonomic constraints in the specification of D. 

Theorem 1. Let cj 7 = £(L)*loq be the pullback of the canonical symplectic form 
on T* D under the Legendre transformation of the prolongation of L: TQ —> M, 
viewed as a map 7: TD —> R. Then T(Yl) viewed as a vector field on TD is a 
Lagrangian vector field for the energy function h: TD — > M. In local coordinates 
(?»£)?)£) forTD, 

(39) h(q,e,q,e) = — q- — e, 
which coincides with Al 



Proof. We carry on the proof in local coordinates. Let (q, e, q, e) be coordinates in 
TD and (q, e, ir,p) coordinates in T*D. The canonical 1-form in T*D is given by 
9q = n dq + pde. The Legendre transform £(7) is given by 

(40) ,-g, 

(41) p= m = W 

then 

(42) ^ = Tq dq+ Tq de 

An straightforward computation shows that 



9 2 7 9 2 7 . 9 2 7 



(43) 



- — —dq k A dq 3 + dq k A dqj + . de k A dq, 

dqkdq-j oq k dqj oe k dq-j 

d 2 j d 2 j 9 2 7 

— de k A dqj + - — ^dq k A dtj + dq k A dej. 



de k ddj dq k dej dq k dlj 

Consider the second order equation (or spray) 



(44) T ^=4q + 4e +C ¥q +D §i 

then 
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(45) 

t{Yl) ^ = ^k {4kdqj ~ ^ dqk) ~ ^h mk - £k^ d£k ~ ^kfi dkk 

d 2 l . j d 2 7 . j d 2 7 ■ ■ d 2 l -j. 

+ q — Tr^qkdtj + - — K^ekdqj - ■= — — Ejdq k - e 3 dq k 

dq k dej oe k dqj aq k de 3 oq k dej 

d 2 j 9 2 7 9 2 7 

+ a . Q . C k dqj + C k dtj + D k dqj. 

dq k dqj oq k dej oe k dqj 

The differential of h is 

ru ( d 2 7 , . d 2 1 d 2 7 , . d 2 7 ■■ \ . 

\dq k dqj aq k dqj oekdqj oe k dqj J 

<9 7 / d 2 ^ d 2 -/ \ <9 7 

+ -^rdqj - dq k + dq k ej - — (%. 

dqj \dq k dej dq k dej J dej 

Equating coefficients of de k in the equation T(Y/ j )lw 7 = —dh one gets 
_d_ (dL\ . _d_ (dL\ _ dL_ 

[ ' dq 3 \dq k ) q3 + dq 3 \dq k ) J ~ dq k 

by using the fact that C k = q k (here the dots means derivative with respect to 
time) wer recover the orignal Lagrange equations 

— ( — \ - — 
( } dt\dq k )~dq k 

Equating coefficients of de k and dq k , sendous identities are obtained. Finally equat- 
ing the coefficients of dq k lead to the linearized equations (R) . □ 

5. Examples 
A few examples will be presented where the above formulation can be applied. 

Geodesic flow. Here the Lagrangian is just the square of the length 

(49) L = -9abq a qb 

where g a b is the first fundamental form of the metric. The linearized equations are 
properly known as Jacobi equations or as equations of geodesic displacement. The 
prolonged Lagrangian is 

, Kf .s dg ab 

(50) 7 = — — e k q a qb + g a bq a eb- 

aq k 

The Lagrange equations give the geodesic equations 

(51) q a + F abc q b q c = 0, 
and the equations of geodesic deviation 
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(52) e a = Rabcdqbqd e c , 

where, T abc is the affine connection or Christoffcl's symbol, Rated is the Riemann 
tensor, and the dots stand in this case for derivatives respect the arc length s |q1|p5|. 

Newton's equations. Here 

(53) L = -m ab q a q b - V(q), 

where m ab is a constant symmetric matrix. The prolonged Lagrangian is 

dV 



(54) 


7 = —5— e a + m ab q a e b 
aq a 


The Lagrange equation 




(55) 


d ( dj\ c>7 
dt \de k J de k ' 


leads to 





dV 
(56) m ak q k = --— , 

oq k 

which are Newton's equations. In the case of the solution being an equilibrium 
point with the configuration a minimum of the potential, these are the equations 
of small oscillations. Lagrange equation 

d f dj \ d-f 



( 57 ) ,, v . r ; •> 

at \oq k J dq k 

then becomes 



(58) makh + K ak e k = 0, 

where 



d 2 V 
(59) K ab = -?JL. 



Plane motions in an arbitrary 2-potential. 

Choose a reference curve and let s = arc length, z — normal distance from the 
curve, and p — radius of curvature; then taking (s, z) as generalized coordinates 
the Lagrangian is || fu], Q 



(60) L=\ 



z 2 + —(p + z) 



-V(z,s). 



P 
The linearized equations obtained by computing the prolonged Lagrangian yields 

•• 3 (s 1 pd 2 V\ 2h 
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Here is worth to say that we choose the "virtual displacements" e s and e 2 to be 
planar ones. The constant of motion h can be expressed as 

■ ■fs 2 dV\ dV 

62 h = sz[— + —)e z + — e s . 

\ p oz J os 

We pinpoint that from equation ( |6l| ) we can see the that a planar orbit is stable, 
against energy preserving perturbations (h = 0), only when the expression between 
parenthesis is positive S |2(| . 

Rotating Lagrangian systems 

Let q — R(t)Q, RBr 1 — il then the Lagrangian of a mechanical system L(q, q) = 
h(q,q) — V(q) transforms into 

(63) L(Q,Q) = ±(Q,Q) + (nQ,Q)-V ef (Q) 

in a rotating frame, where fi is the rotation matrix and where the effective potential 

is 

(64) V ef (Q) = V(Q) + -(1Q, Q), 1 = inertia matrix; 

and (, } is the usual inner product in R™. The prolonged Lagrangian yields the right 
linearized equations: 

(65) Q + 2OQ + W e /(Q) = 

(66) e + 2Qe + Be = 0. 

where B = HessV e f(Q). 
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